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The integration of the storage differential equation a1 present is usually done mostly by 
graphical or numerical procedures. An approach to the analytical integration of that equa- 
tion is the subject of this paper. A new method of fitting the given background curves by 
mathematically tractable expressions is introduced. The storage-outflow discharge relation 
is expressed in the form of a, power function. A general differential equation for water 
storage y' -\-cPy 2 — cy k = is derived, with c and /, constants for the given reservoir, outflow 
shape .'Hid type of flow, and P being the inflow hydrograph. The analytical solutions of this 
equation for P=0, P— constant, and certain P=f(t) are given for the integrable cases (tables 1 
to 3, eqs (12) to (29)). The application of the results obtained is discussed at the end of the 
paper. 



1. Introduction 

The known simple balance equation, with the 
inflow minus tin 4 outflow in a time interval equal to 
the storage change in that interval is 



P-Q= 



(IS 
dt 



(1) 



which is normally derived from the definition itself; 
with P the inflow, Q the outflow, S the storage 
volume, and t the time. 

The two basic differential equations which govern 
the unsteady water movement through lakes and 
channels are the continuity and the dynamic (or 
momentum) equations: 



d(Av) dA 
dx h df = 



=0, 



dv . dv , dh 



(2) 



(3) 



where 



A= cross-sectional flow area, 
fl=mean velocity of flow, 
x= length along the lake or the channel, 
/=time, 

g= earth acceleration, 
A —level of water surface in a cross sec- 
tion with reference to a fixed level, 
and 
S r / =frictional slope. 



*This work was prepared under a National Bureau of Standards contract with 
The American University, with the sponsorship of the Army Map Service. 



The applications of eqs (2) and (3) must satisfy 
the basic assumptions made by De Saint Venant in 
their derivation, and specifically that the unsteady 
flow is gradually varied. 

The eq (3) is omitted in main engineering prob- 
lems, as in the routing of floods through greal reser- 
voirs and along channels of small slopes with un- 
steady (low very gradually varied, so that only eq (2) 
is used. 

Two important problems which have not as yet 
been completely discussed, and are not considered 
here, are: (a) Under what conditions can eqs (2) 
and (3) be used, and (b) when can eq (3) be neg- 
lected without scnsihle error. 

The first term of eq (2), b(Av)/(bx) = (i>Q)f(dx), 
can, for suitable small Ax, be replaced by (AQ)l(Ax); 
for a given time interval At and a reach of length Ax, 
can be replaced by AQ/Ax=(P—Q)/Ax, where P is 
arrival or inflow discharge and Q is the departure or 
outflow discharge in m?/s or in cfs. (Here P and 
Q are used instead of usual / and 0.) The second 
term (dA)/(dt) can be replaced for small ±1 by (AA)/ 
(A0 = (Ail 1 +A^l 2 )/(2A0, where AA 1 =A[ f — A[ and 
AA 2 =A 2 — A' 2 , and A/1, and AA 2 are changes of 
cross-sectional areas at the beginning, and at the end 
of the reach Ax during time At, and A[ and A 2 are 
areas for the beginning, and A[' and A 2 areas for the 
end of time interval A^. Equation (2) can thus be 
expressed as 



p-fr 



_ AAj+AA 2 Ax_ 
2 At = 



AS 
"At' 



where AS is the volume change for the reach Ax 
during time interval At. If now At— >0, this differ- 
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ential equation becomes eq (1), in which S is the 
storage volume for body of water of the considered 
reach Ax, and dS/dt is the rate of change with time of 
this volume. In the case of a reservoir or lake, with 
small velocities along the lake and with inflow com- 
ing from different tributaries, the reach Ax is replaced 
by the water body of the lake. Strictly speaking, 
some end portions of lakes with shallow water and 
high inflow should be considered separately, but their 
influence on the outflow hydrograph can be neglected 
in comparison with the main influence of great 
fluctuating storage space. Equations (2) and (1) 
are identical, and the equation of continuity in its 
form of eq (1) will be called the storage differential 
equation for water in lakes and along channels of 
small slopes with unsteady flow very gradually 
varied with time. 

The analytical integration of the storage differ- 
ential equation would be useful in the solution of 
problems in cases where its terms can be approxi- 
mated by integrable expressions. Some of these 
cases are: (a) Routing of very gradually varied floods 
along large channels of small bottom slope; (b) effect 
of unregulated lakes on flood waves; (c) routing of 
flood peaks through reservoirs for the inflow dis- 
charges higher than that discharge of hydrograph, 
which corresponds to the full capacity of open gates 
and valves; (d) study of the genesis of hydrograph 
and separation of water flow according to its origin 
(surface runoff, underground storage, lakes, channel 
storage, etc.); (e) outflow through partial openings 
in case of rupture of dam with no influence of tail- 
race levels on the outflow hydrograph, when the wave 
movement along the lake created by rupture can be 
neglected; (f) computation of the seepage water out 
of reservoirs, etc. 

Equation (1) serves generally for the computation 
of relations between five functions : 

(1) Inflow hydrograph, P=f 1 (t); 

(2) outflow hydrograph, Q=j 2 {f) ; 

(3) stage hydrograph, H=j z (t) ; 

(4) outflow rating curve, Q=f±(H) ; 

(5) storage function, S=j 5 {H); or area function 

A=j 6 (H), with five variables, Q, P, H, S, t. 

When three of five functions with boundary condi- 
tions are known (three variables can be excluded), 
the eq (1) enables the computation of the relation 
between two remaining variables. 

The usual integration procedures are numerical, 
graphical, combined numerical and graphical, and 
integration by special devices. 

The analytical integration of eq (1) is not usually 
feasible in practice, mostly because of the difficulties 
of fitting easily three known curves out of five by 
tractable and integrable mathematical expressions. 

The most common cases in the application of eq 
(1) have, as given, the following curves: (1) Storage 
or area function (obtained by survey) ; (2) outflow 
rating curve (obtained by gauging, hydraulic com- 
putation, model study, etc.) ; and (3) inflow (or out- 
flow) hydrograph. Two other functions are normally 
to be computed; (1) outflow (or inflow) hydrograph, 
and (2) stage hydrograph. 



The subject of this paper is the analytical integra- 
tion of eq (1) for those cases for which inflow hydro- 
graph, storage function and outflow rating curve can 
be entirely or partly fitted by simple expressions 
which make eq (1) integrable. 

2. Fitting of Mathematical Expressions to 
Given Curves 

2.1. Storage Function 

The storage function, which relates lake or channel 
volume to its level referred to some datum, can be 
approximated either by the function (see references 
[1] l to [5]) of the type 

S=aH m (4) 

or by the polynomial of the type 

S=A +A l H+A 2 H 2 + . . . +A m H m (4a) 

where S= storage volume, H= depth of water above 
a reference level suitably selected, a and m, and 
A , A h A 2 , . . ., ^4 w are coefficients to be determined 
from data of storage function. 

The reference level of eq (4) is normally that of 
zero storage (lowest level of lake or reservoir, or 
river bed, etc.). The level of zero outflow does not 
coincide generally with the level of zero storage, and 
in that case another form of eq (4) is used: 



S=a{H m -H r S) 



(4b) 



where H — difference of the level of zero outflow 
and of the level of zero storage. The coefficients a 
and m in eq (4b) are the same as in eq (4). The 
use of eq (4), as will be shown later, gives a form to 
eq (1) which is integrable in many cases, but eq 
(4b) makes it less simple and more complex for 
analytical treatment. 

Equation (4a) is well-suited to be used as the 
storage function, if enough terms (m+1) are in- 
volved. It must be supposed that for the level 
H=0, the storage is A , so that eq (4a) can be used 
in case the level of zero outflow coincides with the 
level of storage A . If the level of zero outflow is 
changed, the coefficients A , A u A 2 , . . ., A m have 
also to be changed. 

Instead of using eq (4b) repeatedly with various 
values of H to cover the entire range of S versus H, 
a family of eqs (4) may be used with continuously 
changing a and m as function either of level of zero 
outflow or of H C) the distance from level of zero 
outflow to a reference level. The coefficients a and 
m of eq (4) depend on level of zero outflow, or on 
difference of levels H c , and are determined by the 
condition that H=0 of new zero storage coincides 
always with the level of zero outflow. This means 
that the coefficients a and m are to be determined 
for different levels of zero outflow (and therefore also 
of zero storage) as continuous functions. The 

1 Figures in brackets indicate the literature references at the end of this paper. 
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question arises whether eq (4) his the thus conceived 
changing storage function as well as it fits the entire 
function from the lowest to the highest level of the 
body of water. In some eases eq (4), used in the 
manner described above, fits the upper parts of the 
storage curve better than t he entire curve, but in 
other cases the reverse is true. 

Figures 1 and 2 show the second case, and figures 
3 and 4 the first case. Figures 1 and 3 show the 
storage functions S=aH m for the entire range of 
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Figure 1. Storage function S — aH m , with a and m as functions 
of H c , the levels of outflow zero {first example). 

H' 




40 60 80100x10° 



S(m ! 



Figure 2. Storage function of figure 1 in logarithmic scales for 

different values H c of zoo ovtflow, where the storage functions 
from H c up to the highest level are fvtted by the power functions 
S — aH m (first example). 



fl r =depth of water in reservoir at dam, and the 
functions m=ji(H c ) and a=f 2 (H e ) with //, as dis- 
tance of level of outflow zero to the sea level, or to 
the reservoir bottom at dam H' c . Figures 2 and 4 
show the fitting of eq (4) to the parts of storage 
curves for different values H c . When // ( . is ex- 
pressed as H' c , or as distance to the bottom of reser- 
voir, it can be used as dimensionless number h= 
H' c /H m , with iZ m =maximum depth of reservoir in the 
normal operation. 
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Figure 3. Storage function S = aH m , with a and m as functions 
of H c , the levels of outflow zero (second example). 
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Figure 4. Storage function of figure 3 in logai ithmic scales for 
different values H c of level of zero outflow, where the storage 
functions from H c up to the highest level are fitted by the power 
functions S = aH m (second example). 
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The storage functions in figures (1) and (3) are 
given from the reservoir bottom at dam up, and in 
figures 2 and 4 from the level of zero outflow up 
for each of given curves (as straight lines). 

The function m=fi(H c ) decreases with higher 
value of H c , i.e., higher level of zero outflow, and is 
approximately 1 for the highest levels and small 
range of level differences, where the storage function 
is nearly a straight line. For lower H c and for higher 
range of storage levels, the value m is greater. The 
function a=j 2 (H c ) has the opposite characteristics. 

If only one value of H c is to be investigated, a and 
m are fixed coefficients, but if a range of H c is to be 
investigated (dam-breaching problems with various 
breachings and different lowest outflow levels), the 
functions m and a of H c can be easily obtained as in 
the case of figures 1 to 4. 

When the sedimentation of reservoirs has to be 
taken into account, a and m are not only a function of 
level of zero outflow, but also of time (apart from 
sediment characteristics they depend on operation 
practice of reservoir, inflow Irydrographs, etc.). As 
a and m make eq (4) less tractable in case they change 
with time (which is sometimes rather difficult to 
predict), they will be considered here as constants for 
a given reservoir, for given level of zero outflow, and 
for selected time period of reservoir life. 

The coefficient m ranges within the limits 1 to 5, 
but for the majority of natural valleys the range is 
about 1 to 4. When H c is high and the range of 
levels is very small, m is mostly 1.0 to 1.5 and rarely 
greater than 2; for highest range it is 2 to 5. The m 
depends not only on the range of levels, but also on 
the shape of reservoir cross sections. When those 
sections can be fitted by a power function for width 
of type B=2pH s , m is function of s. The value a 
depends on H c , but apart from that, a is higher for 
wider valleys, for lower river slopes, and for lower 
values of exponent s. 

The property of the storage function to be well- 
fitted by an expression of the type of eq (4) can be 
used to compute this function in the case where only 
two points of the relation storage versus level are 
available. These two points can be found in the 
literature for practically any reservoir: capacities at 
highest and at lowest operational levels. Two pairs 
of values (S, H) are sufficient for the computation 
of values a and m in eq (4), thus allowing the deter- 
mination of the curves of figures 1 and 3. Taking 
into account the errors in the survey of reservoirs 
and the fact that the storage function is being con- 
stantly changed by sedimentation and littoral erosion, 
the values a and m determined by eq (4) and the 
figures 1 to 4 are accurate enough for the analytical 
solution of practical problems. 

2.2. Outflow Rating Curve 

The outflow rating curves are relations of the 
departure discharge to the water level in the res- 
ervoir. They can be fitted also in many cases 



either by an expression of type 



(5) 



or by a polynomial, 

Q=b +b 1 H+b 2 H 2 +b 3 FP+. . .+b r H r , (5a) 

where Q= departure discharge, H= level above the 
depth for which Q=b , b and r or b , b u b 2) 6 3 , . . . 
b r are coefficients depending only on boundary 
conditions (outlet shape) and on type of flow (closed 
or free outlets). The eq (5a) can be also conceived 
as a rating function expressed in power series form. 
The range of r in eq (5) for free surface outflow is 
1.5—4.0 (usually 1.5 to 3.0) and for pressure type 
flow about r=0.5 in most cases. The range r= 
0.5—1.5 can be covered also in special cases for 
special shapes of outflows. 

In case two types of flow, viz, flow under pressure 
and free surface flow, occur simultaneously, eq (5) 
takes the form 

Q=b 1 H l »+b& r§ =b 1 (H+H a ) n +b*B», (5b) 

where H 2 =H and H l =H+H a are the hydraulic 
heads (level differences) to be taken for the former 
and latter types of flow respectively, &i and r x 
correspond to closed pressure flow and b 2 and r 2 to 
free surface flow. As r 1 and r 2 are not generally 
integers, eq (5b) is different from eq (5a). 

The use of eqs (5a) and (5b) instead of eq (5) 
makes eq (1) integrable in a much smaller number 
of cases. Equation (5) will be used exclusively 
here. 

2.3. Inflow Hydrograph 

It is practically impossible to fit an entire natural 
inflow hydrograph by a single mathematical expres- 
sion, but it is possible to fit some of their parts by 
tractable and simple functions. 

The simplest forms of function P=f(t) to be used 
for integration are: (a) P=0 (arrival discharge zero 
as dry period of year or as w^ater fully stored in 
upstream reservoirs, etc.); (b) P=P = constant (or 
nearly constant, with low river flow slowly changing, 
regulated constant flow from upstream reservoirs, 
short-term operation of reservoirs, when Row could 
be placed by constant discharge); (c) P=P —ft 
(gradually varied flow); P=P t~ s ; P=P c" /l ; P = 
P t s e~ ft (P , /, and s are different constants in each 
case). 

2.4. Storage-Outflow Discharge Function 

When the lowest outflow level coincides with the 
level of zero storage, or when reference levels for 
measuring H in both eqs (4) and (5) are the same 
(adjustment to this condition is always possible by 
taking the corresponding coefficients a and m from 
their functions, fig. 1 to 4), the elimination of H 
gives for the storage-outflow discharge function, or 
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relation $=j'(Q), 



S= 



h- 



Q»/'=± Q" t 



(6) 



where c = b {m/r) la and n=m/r. 

The coefficient n for free surface flow is in the 
range }i to 4, and for closed pressure flow 2 to 10, 
or a little greater. 

The use of pairs of eqs (4) and (5a), (4) and (5b), 
(4a) and (5a), (4a) and (5b), (4b) and (5), (4b) and 
(5a), (4b) and (5b) gives a more complex relation 
S=f(Q), in most cases in parametric form with II 
as (he parameter, than the use of (4) and (5) or (4a) 
and (5). However, in this last case 

S=a +y Q^+~ <2>/'+ . . . +1 Or*, (6a) 

c l ^2 Cm 

which makes eq (1) less tractable than the use of 
eqs (4) and (5). 

As eqs (4) and (5) fit well many storage and out- 
flow functions, for different value of II c , and as 
eq (6) allows the integration of eq (1) in many cases, 
eqs (4) and (5), and (6) will be used in this study. 

3. General Type of Differential Equations 

With the introduction of eq (6) into eq (1), one 

obtains n Q n ~ l dQ+c(Q-P) dt=0. 

The substitution y = Q~ n and the replacement 
I r 

k=--(2n—l)/?i = 2 =2 — - (hen gives the following 

expression : 



y'-{-cPy 2 —cy k =0. 



(7) 



This is the general differential equation for storage 
reservoirs and storage channels with no artificial 
control of outlet flow, under conditions discussed 
previously in this paper. For free surface outflows 
the range of k is usually from about —2 for n=}{ to 
fc=% for ?i = 4, and for the closed type outflows the 
range of k is from % to 1 %o or a little greater, but 
less than the limit value k=2. Practically, the range 
of k is from —2 to +2, where k=+2 is not possible. 

The relations of k or n to m and r are given in 
figure 5. The most common range for both m 
(1 to 4) and r (0.5 and 1.5 to 3) is specially empha- 
sized for reservoirs in river valleys. In most 
practical cases &>0, and usually &>l/2. 

The eq (7) could also be expressed as 

y'+cy 2 (P-l)-ey r/m = Q. (8) 

The general form of eqs (7) and (8) is 



y f + Mi/ l +Ny g =0 



(9) 



where M and N are functions of t; h and g are any 
real numbers. By the substitution y=Z {lnh ~ l)) the 
general form is reduced to 



h+g-2 



Z' + (h-l)MZ 2 +(h-l)NZ "- 1 =0 



(10) 




6.0 



FiGt re 5. \ r alues n of storage-outflow discharge relationship 

S = —Q", and values fc=2 , in .elation to m and r (expo- 

c ' n 

nents of storage and of rating curve, expressed in power form). 



which is eq (7) under conditions that h— l=c, 
N=-l, M=P, and (h+g-2)/(h-l)=k. 

The general analytical solution of \\u x ordinary 
differential equations of the type (7) or (10) for any 
/• and for any expression of P does not exist. 
Chiellini [6] has shown that the solution is possible 
for fc=3 if the expression (M'N—MN f ) is different 
from A/ 3 by a const an! factor. When k is a positive 
and integer number higher than 3, Luigi Conte 
[6, 7] has given the solution for special cases. As k 
is always less than two here, their procedures can- 
not be used for eq (7). 

Both facts, (a) that k is a number derived, from 
continuously changing natural conditions and is 
rarely an integer, and (b) that the function P=f(t) 
is not simple and must be (it ted by integrable 
expressions, make the analytical integration of eq 
(7) possible only in special cases. For k —2, —1, 
0, and 1 and some rational numbers the integral ion 
is simpler for some types of P, but the integration 
for intermediate values of k and for the majority 
of P functions has to be performed numerically or 
graphically. The analytical integration for special 
values of k and for selected functions P=f(t) is 
useful, because the resulting functions give the type 
of curves which can be expected for intermediate 
values of k. Or, when k and outflow functions are 
known, some conclusions can be made about the 
inflow functions. The type of functions obtained 
in integrable cases can serve for control of curves 
obtained by other procedures. 
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4. Integration in Special Cases 

Equation (7) will be analyzed for some values of 
k and some types of the function P. The values 
of k taken will be — 2, — 1, 0, 1 and, tentatively, 

— 3/2, —1/2, 1/2, 3/2 along with some other rational 
numbers. 

4.1. P=0 



Equation (7) becomes 

2/'-q/*=0, 



(11) 



with y=yo=Qo~ n for t=0, where Q is outflow dis- 
charge at the time t=0. The solution of eq (11) 
with y=Q~ n is then 



e— i =«s- 1 — (i— |) c*. 



(12) 



This is the general expression for outflow hydro- 
graphs from storage reservoir with no inflow into the 
reservoir. The solution can be obtained from eq 
(12) for each value of n except for n=l (k=l), 
but eq (11) gives solution in this case directly 
(table 1). Therefore, for P=0, or when the constant 
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inflow discharge is small and can be neglected in 
comparison to the outflow discharges, the analytical 
integration of eq (7) gives solutions for any k(n), 
rational or irrational. 

Equation (12) can be written in the form: 

Q M l - n -^<f) A (13) 

This shows that the decrease of discharge ratio 



Q/Qo depends on both t and Q , except for n=k=l 
(table 1), in which case it depends only on t. 

The analytical expressions of eq (12) for some k 
(or n) values are given in table 1. Figure 6 gives 
curves for n=\ ) 1, f, 2 and 3, for given values of 
a and b, with c—b n ja. The values a of eq (4) and b 
of eq (5) are taken constant for all five curves (for 
easier computation in examples), though they change 
when m and r (and so n) changes. For n=2 the 
outflow hydrograph is a straight line, which corre- 
sponds, among other cases, to the cylindric reservoir 
(m=l) and to the outfhyw under pressure (r=i). 
The same type of linear hydrograph is obtained for 



m = 3 and 



which is the case for some deep 



lakes with outflow function Q=bH%. For river 
outlets from lakes, with r>| and ra<3, n<2, all 
hydrographs are concave upwards. For n=l, the 
solution of eq (11) is an exponential function. All 
functions for n<\ are asymptotic with the time axis. 
Hydrographs for n>\ do not have this asymptote, 
but cut the t— axis for finite values (for n = 2, 
t = 2Q /c). For the range 2^>n^>l, the hydrographs 
are power functions and concave upward, but have 
finite values for £=0 (are not asymptotic to the t 
axis). For n<0 the curves are hyperbolic functions, 
asymptotic and concave upward. 

All outflow hydrographs for very small water levels 
above the bottom of outflow orifice become asymp- 
totic, because their outflow functions change due to 
the change of type of flow or to some secondary 
effects. In closed-type outflow the function changes 
from r=i to about r=f as soon as the water 
level drops to near the upper edge of the outflow 
orifice, and changes further for smaller outflow 
heights due to surface tension and other secondary 
effects. The exponent r increases so that n=m/r< 1. 

The shapes of outflow hydrographs depend thus 
on the ratio of exponents of storage and outflow 
functions, and they can be: hyperbolic (7i<l), 
exponential (n=l), power function with exponents 
higher than 1 (2>n>l), straight line (n=2) and 
power function with exponents lower than 1 (n^>2). 
They show that the natural conditions can produce 
different outflow hydrographs on the same river, 
depending on the lake and channel characteristics 
and outflow shapes. 

The free spillways of reservoirs, or free outflows 
of shallow lakes, where the range of storage fluctu- 
ation is small, m=1.0 to 1.3, and r>f, have n<l, 
and all outflow hydrographs are hyperbolic asymp- 
totic functions (fig. 6). For the higher range of 
levels of outflow, m is greater, and the outflow 
hydrographs are power functions whose graphs are 
convex upward or downward. The smallest value of 

n for closed-type outflow is 71=;— =2, so that all 

hydrographs are power functions convex upward. 

In case eq (6a) is used instead of eq (6), the follow- 
ing general solution of eq (7) for P=0 is obtained: 



*«s 



m 



i CJm—r) 



(m—r m — r \ 

Q7~ -Q~ i 



(14) 



48 
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Figure 6. Outflow hydrographs in their recession parts from 
Q to zero, related to different values k and n, in the case where 
inflow is zero. 

which for the case m> 3 is more complicated than 
eqs (12) or (13), particularly since eq (14) gives t 
as a function of Q rather than the more convenient 
inverse. 

If both types of flow occur simultaneously, eq (5b) 
has analytical solutions in a much more limited 
number of cases. 

4.2. P=P =Constant 

In the case of constant or nearly constant inflow 
into the reservoir, with P=P , eq (7) becomes 



y' + cP y 2 — cy k =Q, 
and the separation of variables gives 
dy 



J 



P^-y* 



-ct=K 



(15) 



(16) 



with K= integration constant. 

The indicated integration of eq (16) cannot be 
carried out analytically for arbitrary k, contrary to 
the case of eq (11). Equation (16) has, however, 
analytical solutions for many rational numbers, which 
is useful, since it is possible to approximate k by a 
rational number of type gjh, where g and h are 
integers. Equation (16) can be written in the form 
of the binomial integral 



/ 



y- 2 (I > o-y k '- 2 )- , <ly+ct=K. 



(17) 



The binomial y 2 (P„—if 



is of the form 



y s (ay r +b)» 



(18) 



which can be reduced to an algebraic function in case 
any ono of the three Qumbers 



P, 



*±1, 
r 



and 



8+1 



-v 



is an integer. As p = -1 is an integer, when k is a 
rational number, eq (16) can always be reduced to 
an algebraic function. 

The most general solutions of the binomial integral, 
if the obtained algebraic functions are integrable, 
are in the form of sum of logarithmic, algebraic and 
arctan terms. In some cases the a ret an term does 
not appear in the solution, as is the case also with the 
algebraic term, but the logarithmic term is always 
present. The solutions of eq (16) will be analyzed 
here for some values of k. 

By substitution of a new variable u in eq (16), 



u=y 1 < h -P Q 2h -e, 



(19) 



the following general binomial expression of eq (16) 
for k=g/h is obtained: 



g—h 



with h — g 
For 



hP 2h -o J u h - g - 1 (u 2h - g —l)- l (!u+ct = K (20) 
1 and 2h — (j integers. 



f(u)=ju h -'-* 



i„ 



2h-g_ 



\)-\lu 



(21) 



the constant Kfort=Q (u=u ) isK=hP (g - h)/m - g) X 
/0/ u ), and the general solution has the form 



t=Rlf(u )-f(u)] 



(22) 



h 



withR=-P Q (g - h)/m - g) . 



The following values and ranges of q and h are 
analvzed: (1) g=2h) (2) h<g<2h; (3) g<h\ (4) 

(1) g=2h. For k=2, eq (20) lias no practical 
significance. 

(2) h<g<^2h. The value s = 2h—g is always a 
positive integer number, and j=h—g—l is always a 
negative integer number or zero. Equation (20) is 
integrable if s=2h— g<4 or .s = (i, when the rational 
function can be separated in integrable terms, what- 
ever the value of j=h — g — 1. The solution is 



t~- 



h 



r 



Ws) -1 



lf(u ) -/(«)]. 



(22a) 
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(a)s=2h—g=l, n=h. 



/(tO=ln 



u ^u r 2u 1 ^ 



(21a) 



1 (n— l)^" 1 

with u = P Q- 1 . 

This case covers the following values of n and k: 

n: 1, 2, 3, 4, 5, 6, 7, etc. 
k: 1, 3/2, 5/3, 7/4, 9/5, 11/6, 13/7, etc. 

Some of these cases are given in table 2. 
For n=l (k=l) only the logarithmic part of eq 
(21a) is present and eq (22a) gives: 

Q=Po+(Qo-Po)e- ct (21b) 

(b) s = 2h—g=2, n = h/2. 

,, . . u-1 ,2,2, * 



2(n—l)u 2(w - 1) 

(21c) 

with u=P 1/2 Q~~ 1/2 . This case covers the values: 

7i : 3/2, 5/2, 7/2, 9/2, etc. 
fc: 4/3, 8/5, 12/7, 16/9, etc. 

Some of these are given in table 2. 

Table 2. P= P = Constant ; t=R [f(u Q ) —f(u)] 



n 


?? 


k 


P 


w 


f(u) 


1 

4 


0. 25 


-2 


c 


plQ-4 




u—1 , „ 
In — — -4-2 arctan u, 
u-\-\ 


1 

3 


0.33 


-1 


1 — 2 
- P 3 

C 


pi Q - 3 




1 (u-iy /2w+l\ 
In -— — f- +4 arctan ( — — ) 


1_ 
2 


0. 50 





c 


P2 Q-§ 




1 w- 1 


2 
3 


0.67 


1 
2 


2 -1 

-P 3 

c 


P3 Q~J 




in-^i -4 arctan (*±p) 

V?+U+l V 3 / 


3 

4 


0.75 


2 
3 


- 3 p-i 

c 


P*Q~l 




In — — —2 arctan u 

u-\-\ 


1 


1.00 


1 


1 

c 


PoQ" 1 


In — -; ^=Po+(<?o-Po)c-«« 


5 
4 


1.25 


6 

5 


M 

c 





?/4-l 4 

In H 1-2 arctan u 

u—1 u 


4 
3 


1.33 


5 
4 


4 1 

-P3 

c c 


P3 Q~3 




1 . (m-1) 2 . 3 . . . /2u + l\ 
2 1 %.*+i<+l + 4 +2wCtan ( 3 ) 


3 
2 


1.50 


4 
3 


c 


P2 Q-l 




1 w ~ ] , 2 


.5 
3 


1.67 


5" 


-'pi 

c 


P3 Q~3 




WVL+?* 2arctan( 2M + 1 ) 


2 


2. 00 


3 
2 


c 


PoQ-> 


ln — 

u 11 


3 


2.33 


11 


7 4 
-/>3 
c 


1 _ 1 

P3 Q 3 




in <" S - 1)8 +' 


5 

2 


2.50 


8 
5 


C 


1 _ 1 
P2 Q 2 




1 M_1 + 2 4- 2 


3 


3.00 


5 
"3 


-V 

c ° 


Po<?- J 




7 
2 


3.50 


12 


7 5 
-P2 
c 


] _ 1 

pi Q 2 




. m-1 , 2 , 2,2 

ln ?/ +i + /7 + 3 ? 73 + 5¥J 


4 


4.00 


7 
4 


4 3 

c 


PoQ- 1 


ln ~ir+u+2^w 



(c) s = 2h—g=3, 7i=h/3. The solutions for three 
values of h and 71 (A=4, 7?,=4/3; A=5, 77=5/3; 
A=7, 7i==7/3) are given in table 2. 

(d) s = 2h — g=4, n=h/4t. For h=5, n = 5/4 the 
solution is given in table 2. 

(3) g<Ch. The analytical integration can be 
carried out when — 4<2A — g=s<4, or #/2 — 2<h 
<#/2 + 2. For h>g/2 + 2 and h<g/2-2 the ex- 
ponent s=2h— g in eq (20) is higher than 4, so that 
the algebraic function u g ~ h+1 (u 2h ~ g — 1) can not be 
separated into integrable components except in 
special cases. The solutions for n . = 3/4 (fc=2/3), 
?i = 2/3 (£=1/2), n=l/3 (^=-1) and 71= 1/4 
(£=—2) are given in table 2. 

(4) gf=0, A=l; n=l/2 (k=0). The solution is 
given in table 2. Figure 7 gives five curves for the 
values 71 : 1/2, 1, 3/2, 2 and 3, as in figure 6. All five 
are asymptotic to the value P =0-20, but for t nearer 
to zero they approach the curves of figure 6. 



4.3. P=f(t) 

The solution of eq (3) for four types of P=i(t) 
can be obtained for special values of k, but the 
general analytical solution such as was obtained in 
the case of eq (13), or in many special cases of eq 
(22), is not possible. 
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Figure 7. Outflow hydrographs in their recession parts from 
Q to inflow P=P Q = constant, related to different values k 
and n. 
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a. k-l (n = l) 

Equation (7) becomes 



(23) 



which is a Bernoulli equation, and lias the general 
solution for y=Q~ 1 : 

Q=e~ ct (K+c fPe ct dt\ (24) 

The solutions for four types of function P=f(t) 
are given in (able 3. 



Table 3. P = f(t) 
k=l, n=l, s= integer 



p=f(t) 


Q=h(t) 


r= p Q -ft 


(r.-f t+ i) + (,,-n-^ 


p=pj-° 


[«?.-t(-l)-.>]e-+t- ! (- 1 ){l-^f-...+(-l)^] 


P=P e-f* 


(-£)-+(a)- 


Pn= r 


[*"<-nS^]— ^cSy- 




[l-l!(r-/)/+2!(c-/)-'/- > - . . . +.v!(r -/>7*|r '' 



Such analytical solutions can be obtained under 
condition thai I lie integration of 

| Pe ct dt (25) 

in eq (23) can be carried out in closed form. 

b. * o ( n =l) 

{Munition (7) becomes 

y'+cPy 2 -c=Q (20) 

which is a special type of the Riccati equation, with 
y=Q~ 1/2 . As two particular integrals of eq (26) are 
known, viz, y = Q ~ l/2 for /=0, and Q=0 or Q=P 
for t infinite, the solution of eq (2(5) is possible. 
The substitution Z=y—y gives 

Z' +cPZ 2 +2cPy Z+ ( y ' +cPyl-c)=0 

and eq (26) becomes the Bernoulli equation 

Z'+cPZ 2 +2cPy Z=0 (27) 

with t he solution 



'"J y 



1 2cy Pdt 

--e J (K+c 



r -2cy I Pdt 

\Pe J d 



t) (28) 



y— yo 

The constant K is to be determined from the known 



value of Q for / infinite. 

The solution for P=P e~ ft is 



where 



/ ,,-//_ 1 \2 

g=2cP /fy/Qo. 



(29) 



The analytical integration is not possible for tin; 
three other P functions of table 3, because the integral 
fe atS dt can be obtained, in general, only in power 
series form. 

This analysis shows that the analytical integration 
for simple expressions P /(/) can be done for k=l 
and for fc=0 in some cases. 

5. Discussion 

The analysis of the feasibility of lit ting the given 
curves, as parts of differential equation for water 
storage, by mathematical expressions and its inte- 
gration by analytical procedures in special cases, 
show some characteristics of outflow hydrographs. 

The shape of outflow hydrographs from Lakes and 
channel reaches is highly influenced by the ratio of 
exponents of storage function and outflow rating 
curve, both of the power type y = sx p . 

The analytical integration of the storage differen- 
tial equation can be performed for any ratio of ex- 
ponents in case of zero inflow, and for many ratios 
of exponents in case of constant inflow. The inte- 
gration is, however, limited to a small number of 
ratios and inflow functions, when the inflow changes 
with time. 

The usual procedure, to plot the recession part of 
river-flow hydrographs on semi-log paper and to fit 
it by a straight line, shows that in the majority of 
cases the fitting of a straight line is a rough approxi- 
mation. The results of above analytical integration 
show that the straight line can be used for accurate 
fitting only if the ratio of exponents of storage func- 
tion and outflow^ rating curve has a mean value near 
unity (n=m/r=l) and for zero or constant inflow. 
This can be considered a special case. 

Though the discharge hydrographs of underground 
water in a river basin can be (theoretically) fitted in 
many cases during the recession period of no water 
supply to the underground by an exponential func- 
tion of type Q e~ ft , the water storage either in lakes 
or along river channels influences highly the shape 
of river hydrographs, so that it becomes less and less 
of the pure exponential type. The 4 higher the ratio 
between the effective storage (storage which influ- 
ences the outflow hydrograph) of lakes and channels 
to the total river flow during the period of recession 
flow, and the higher the departure of the ratio of 
two exponents from unity, the greater is their influ- 
ence on the shape of river hydrographs and the Larger 
is the departure from the original hydrograph of the 
recession flow from the underground. 

As the ratio n = m,/r changes along the river 
channel, the accurate flood routing procedures must 
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take that fact into consideration. 

The described analytical procedure is useful for the 
study of outflow hydrographs in case of rapid 
openings in water bodies (breaches of dams, rapid 
openings of gates and valves, breaches of channel 
walls or levees, etc.). The assumption of simple 
shape of openings and of simple inflow hydrographs 
(zero, constant inflow) is usual in case of dam 
breaches. As the study of dam breaches has to be 
made for many openings of quite different dimensions 
and shapes for different inflow discharges, the 
analytical integration as given above can have some 
advantages in the computation of hydrographs in 
comparison to the standard graphical or numerical 
procedures of integration. 

The difficulties in fitting the background curves 
by tractable mathematical expressions and the 
difficulties of integrating analytically the storage 
differential equation, and especially in case the inflow 
changes with time, explain why the usual graphical 
and numerical procedures have taken so common a 
place in the engineering practice. Nevertheless, the 
analytical treatment of the storage differential 
equation can be useful for some problems, can serve 
for better understanding of relations, and can save 
efforts in problems, as, for instance, in case of outflow 
hydrographs out of dam breaches. 



The assumptions usually taken for some breaches 
and for inflow discharge, and the accuracy of basic 
data in the routing of waves created downstream, 
justify the use of the procedure given for the fitting 
of mathematical expressions to the basic and given 
curves and the use also of the analytical integration 
of the storage differential equation for the com- 
putation of outflow hydrographs. 
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